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METHOD FOR NEAR OPTIMAL JOINT CHANNEL ESTIMATION 
AND DATA DETECTION FOR COFDM SYSTEMS 



Related Applications 

[001] This application claims the benefit of priority of Provisional Application No. 

60/241,195, filed on October 19, 2000, and of Provisional Application No. 60/180,799, filed 
on February 7, 2000, and is related by subject matter to U.S. Patent 6,477,210, entitled 
"System for Near Optimal Joint Channel Estimation and Data Detection for COFDM 
Systems," by the inventors of the present application filed concurrently herev^ith. 

Field of Invention 

[002] The invention relates generally to communications and particularly to a method and 

apparatus for near optimal joint channel estimation and data detection to improve channel 
tracking and, thus, improve link robustness. 



Background of the Invention 
[003] The rapid grovrth in the use of the Internet and the increasing interest in portable 

computing devices have triggered the desire for high-speed wireless data services. One of 
the more promising candidates for achieving high data rate transmission in a mobile 
environment is Orthogonal Frequency Division Muhiplexing (OFDM), v^hich divides the 
wide signal bandwidth into many narrow-band subchannels, which are transmitted in 
parallel. Each subchannel is typically chosen narrow enough to eliminate the effects of delay 
spread. Coded OFDM (COFDM) systems, which combine both OFDM and channel coding 
techniques, are able to improve the performance further by taking advantage of frequency 
diversity of the channel. 

[004] Though both differential and coherent demodulation can be applied in a COFDM 

system, the latter leads to a performance gain of 3 to 4 dB in signal-to-noise ratio (SNR) with 
accurate chaimel estimation. Chaimel estimation techniques realized by a frequency-domain 
filter using Fast Fourier Transform (FFT), followed by time-domain filters for a COFDM 
system with Reed-Solomon (RS) coding have been proposed. These channel estimation 
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techniques, while good, did not provide the near optimal channel estimation required for 
data-decoding with improved channel tracking capability for reliable link performance even 
under high user mobility and/or high RF carrier frequency. 

Summary of the Invention 
[005] An advance in the art is achieved by realizing that an iterative approach for 

developing a channel estimation which employs current information, e.g. decoded data and 
newly received signal, yields a better charmel estimation. At each time epoch n , the 
improved receiver does not proceed to the time epoch n + 1 immediately. Instead, the 
receiver uses the newly generated channel information for the next time epoch, n+1, which is 
derived based on the current information, to substitute for the current channel estimate and 
repeats the calculation for the current channel estimation and the decoding of the transmitted 
symbols. In other words, the same procedures are performed twice. Thus, an OFDM 
receiver of the present invention estimates the channel by receiving a multicarrier signal, 
transforming the received signals by applying Fast Fourier transformations, and estimates the 
channel characteristics of the multicarrier channel using iterative processing, and in the 
process also decoding the transform multicarrier signal. 

Brief Description of the Drawings 
[006] FIG. 1 A is the physical configuration of a system for transmitting and receiving 



[007] 
[008] 
[009] 



signals; 

FIG. IB is a block diagram for a baseband receiver at time instant n; 

FIG. IC is a flowchart for a related sub-optimal channel estimation method; 

FIG. 2 is a comparison of RS and convolutional codes, where the maximum Doppler 



[010] 



frequency is 40 Hz;. 

FIG. 3 is a comparison of RS and convolutional codes, where the maximum Doppler 



[Oil] 



frequency is 200 Hz; 

FIG. 4 depicts the performance of the K=9 CC at different maximum Doppler 



frequencies; 
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[012] FIG. 5 depicts the performance of the iterative approach at 200 Hz maximum Doppler 

frequencies, where K=9; 
[013] FIG. 6A is a flowchart for the iterative processing; 

[014] FIG. 6B is a flowchart for the iterative backward processing; 

[015] FIG. 7 shows the performance of iterative backward-processing approach at 200 Hz 

maximum Doppler frequency, where K=9; 
[016] FIG. 8 shows the performance of iterative backward-processing approach at different 

maximum Doppler frequencies, where K=9; 
[017] FIG. 9 is a comparison of systems having different time-domain FIR estimators, 

where K=9; 

[018] FIG. 10 shows the average retransmission probability as a function of occupancy per 

sector. Eight slots and three RF carriers (< 2.5 MHz) are reused in every base station, each 

with three sectors, using Dynamic Packet Assignment (DP A). In the original (non-iterative) 

method, K=5, 40 Hz Doppler or K=9, 125 Hz Doppler; 
[019] FIG. 1 1 shows the average delay of delivered packets as a function of throughput per 

base station, each with three sectors. Eight slots and three RF carriers (< 2.5 MHz) are 

reused in everywhere using DP A. In the original (non-iterative) method, K=5, 40 Hz 

Doppler or K=9, 125 Hz Doppler; 
[020] FIG. 12 depicts the average retransmission probability as a function of occupancy per 

sector. Eight slots and three RF carriers (< 2.5 MHz) are reused in every base station, each 

with three sectors, using DP A, where K=9, 200 Hz Doppler; and 
[021] FIG. 13 shows the average delay of delivered packets as a function of throughput per 

base station, each with three sectors. Eight slots and three RF carriers (< 2.5 MHz) are 

reused in everywhere using DPA, where K=9, 200 Hz Doppler. 



Detailed Description of the Preferred Embodiments 
[022] The physical layer configuration of the system for near optimal joint channel 

estimation and data detection for COFDM systems is depicted in FIG. 1 A. An exemplary 
transmitter 105 is shovm on top and an exemplary receiver 140 is shovm on the bottom of 
FIG. lA. At transmitter 105, a data stream is accepted by a convolutional encoder 1 10, 
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which convolutionally encodes the data stream. The encoded data stream may then 
optionally be forwarded to an interleaver 1 15 for interleaving. If the encoded data stream is 
interleaved in the transmitter 105 then the receiver 140 must correspondingly de-interleave 
the encoded data stream. After interleaving, a modulator 120, for example a QPSK 
modulator 120, modulates the encoded (and optionally interleaved) data stream, which is 
then forwarded to an inverse Fast Fourier transformer 125, to subject the modulated encoded 
(and optionally interleaved) data stream to inverse Fast Fourier transformation. The 
transformed modulated (and optionally interleaved) encoded data stream (signal) is then 
transmitted, in the present invention, over the air via RF unit 130 and antenna 135. 

[023] Correspondingly, receiver 140 accepts multicarrier transmitted signals (data streams) 

via antennas 145 and RF units 150 and subjects the received multicarrier signals to Fast 
Fourier transformation using Fast Fourier transformers 155. These transformed signals are 
concurrently fed into channel estimator 165 and demodulators 160, for example QPSK 
demodulators. The demodulated transformed signals are combined in maximum ratio 
combiner 170. The combined demodulated transformed are then optionally de interleaved 
using de-interleaver 175. The combined demodulated transformed (and optionally de 
interieaved) signal is then decoded using Viterbi decoder 180. The decoded combined 
demodulated transformed (and optionally de-interleaved) signal is then fed back into channel 
estimator 165, which forwards chaimel estimations, which are added to the transformed 
signals that are forwarded to demodulators 160. 

[024] FIG. IB shows the baseband processing, in particular, the iterative nature of the 

receiver portion of the system disclosed herein for near optimal joint charmel estimation and 
data detection for COFDM systems. Chaimel estimator 165 accepts transformed signal 190. 
Chaimel estimations 198 are fed back into channel estimator 165. Channel estimations 194 
are fed into decoder 185, which comprises maximum ratio combiner 170 (shown in FIG. 
1 A), optional de-interleaver 175 (shown in FIG. 1 A) and Viterbi decoder 180 (shown in FIG. 
1 A). Channel estimations 195 are fed into decoder 185 via demodulator 160 (shown in FIG. 
1 A, but not shown in FIG. IB for clarity and to highlight the iterative nature of the system), 
which demodulates the transformed signal using channel characteristics. Decoder 195 also 
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accepts transformed signal 190. Decoder 185 outputs signal 192, which is fed back into 
channel estimator 165. 

[025] An OFDM signal is divided into a number of subchannels. By way of example, an 

OFDM signal bandwidth is divided into 120 6.25-kHz subchannels with QPSK modulation 
on each subchannel. At the receiver, the demodulated signals from two receiving branches 
are combined using maximal ratio combining and then decoded. With a symbol period of 
200 s (including a 40- // s guard interval) and 1/2-rate coding, a maximum information rate 
of 600 kbps can be achieved in a 750-kHz bandwidth (about 800 kHz including guard bands). 
The information rate is calculated by dividing the 120 subchannels (tones) by the 200 /i s 

period to obtain 600 kbps. 
[026] For purposes of example for the present invention, 1/2-rate convolutional codes (CC) 

are considered. The results with 1/2-rate Reed-Solomon (RS) code based on Galois-Field 
(64) (GF(64)) are compared. The size of a code word is the same as that of an OFDM block 
(an OFDM symbol of 200 s and 120 subchannels). To achieve coding gain with inherent 
frequency diversity in OFDM, a simple interleaving scheme is applied. For both RS and CC 
cases, the first 120 bits of a code word are assigned to the in-phase component and the rest to 
the quadrature component. To gain additional randomness within a code word for the CC 
case, each 120-bit group is interleaved over subchannels by an 1 1-by-l 1 block interleaver 
(without the last bit). 

[027] In the simulations, the wireless channel, as a Rayleigh-fading charmel, with a two-ray 

muhipath delay profile is modeled. Good performance for impulse separation as high as 40- 
ju s can be achieved; a. 5- /us impulse separation in the numerical results is considered. 

[028] For the performance with respect to channel variations, maximum Doppler frequency 

up to 200 Hz, which is reasonable for most vehicular speeds, for a possible RF carrier 
frequency around 2 GHz is considered. To demonstrate the advantage of the proposed joint 
detection methods, results at a maximum Doppler frequency as high as 500 Hz corresponding 
to a scenario in which the wireless system uses a higher carrier frequency, e.g. 5 GHz are 
presented. 



5 



Chuang 2000-0071 

[029] In the medium access control (MAC) layer, a frequency reuse is considered with 

dynamic resource management, e.g., Dynamic Packet Assignment (DP A), to achieve high 

spectral efficiency for packet data access. 
[0301 A simple analysis to highlight the ideal or optimal joint charmel estimation and 

maximum likelihood (ML) decoding scheme indicated in FIG. 1 A for the case of M = 2 

receiving antennas is now presented. 
[031] At a diversity receiver, the signal from the m th antenna at the k th subchannel and 

the n th block can be expressed as 

where a„ , Knk >v„ „ ^ are the transmit signal, channel response and additive Gaussian 
noise, respectively. 

[032] For convolutional codes, because the size of a code word is the same as that of the 

OFDM block, (1) can be rewritten as 

V«=H;„,«^„+w,„, (2) 

where, if there are subchannels, H„„ =J/ag(A„ 2 v..,^„,„,jc,) ^ c„ is the transmitted 

code word at time epoch n , and the rest of the vectors are similarly defined. 
[033] Assume that the number of code words is A/' , we introduce the following notations, 

C — [C| J , 

H;„ = t//ag(H^ ,,H^ 2,...,H„ ^) , (3) 

_ r r T r-if 

[034] At the receiver, the objective is to solve a maximum likelihood (ML) problem 

c = argmin[minXlk-H.c||^], (4) 

with a constraint on channel response 

AHJ = 0, (5) 
where LQ is a constraint function. In a wireless environment, this constraint can be 
simplified to be 
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X;B„,d(H..„.;) = 0, (6) 

where the length of the channel memory is OFDM symbol durations, B„ / are 
coefficients determined by the correlation between channel responses at the time epochs n 
and n-1, which is a function of the Doppler spectrum of the channel, and dQ is a vector 

function defined by d(H^ „ ) = „ , , h^^^^ h^^ ^.^ f . 
[035] The optimal solution of this ML problem can be obtained by exhaustive search. It 

requires solving the mean square error (MSE) 

MSE(c) = minXK-H„c||\ (7) 

for any possible c with the channel constraint (6). Then, 

c = arg min MSE(c) . (8) 

c 

[036] After obtaining MSE(c) , the corresponding channel estimate H^(c) can be found. 

Consequently, the optimal approach for estimating channel response requires the knowledge 
of the entire set of x and c . 

[037] Another observation from this ML receiver is that the channel estimation results 

is not a direct output of the detection process and hence, channel estimation which calculates 

explicitly may not be necessary in theory. However, for other required parameter 
estimation, such as timing and fi:equency synchronization, a known data sequence is usually 
transmitted in the beginning of a group of OFDM blocks. This knovm data sequence, also 
called a synch word or a unique word, can be used as a training sequence in (7) to obtain 
initial channel estimate explicitly without resorting to blind detection. This initial channel 
characteristic is helpful for solving this ML problem with better numerical stability and 
tracking property. This initial channel estimation can be easily solved in the frequency 
domain by first taking FFT as shown in FIG. 1 A. 

[038] One related method and system is Serial Number 09/089,862, METHOD AND 

APPARATUS FOR CHANNEL ESTIMATION FOR MULTICARRIER SYSTEMS, which 
was filed June 3, 1998, and is commonly held and incorporated herein by reference. The 
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near optimal joint channel estimation and data detection method and system of the present 
invention was bom from the research that resulted in that application and the subsequent 
determination that improvements could be made in the channel estimation. 

[039] The sub-optimal approach of the related system and method is now outlined. Because 

of the formidable complexity of the optimal ML receiver, some sub-optimal solutions are 
widely used in practice. The related sub-optimal solution is to divide the ML problem into 
two parts, channel estimation and coherent decoding. Then, the problem can be solved by 
iteratively estimating channel and decoding in the forward direction (in time). 

[040] At a time instant n , given a channel estimate H„ „ , initially obtained by using training 

sequence (in the frequency domain dividing the transfer function of the received signal by the 
transfer function of the known data), the maximum likelihood (ML) problem 
c„ = arg min ^\K.n - ^m^'^n f (9) 

[043] can be solved. Then, the reference for channel estimation is 

H„.„ = argmin^||x„ „ - H„ „c„ ||' . (iq) 

[046] Finally, the estimate for the time instant « + 1 is obtained by solving a linear constrained 

equation (6). Considering a stationary channel with fixed maximum Doppler fi-equency, the 
coefficients B„ ^ are independent of n , and can be written as , which are used as the 
coefficients in the FIR filter to track channel variations. Consequently, a simplification of (6) 
with only previous references for prediction-type estimation leads to 
Ml 

lB,d(H„,.,.,)-d(H„.„.,) = 0, (11) 

[049] where is the number of taps of an FIR filter and are preset coefficients designed to 
achieve the minimum mean square error (MMSE) of estimation. This MMSE estimator can 
be realized by a frequency-domain filter using the Fast Fourier Transform (FFT), followed by 
an -tap time-domain filter, 

B,=6,F%F, (12) 
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[052] where is the time-domain filter coefficient, F is the FFT matrix, B^^ is a diagonal matrix, 

and F"*B^F is the frequency-domain filter. 
[053] The MMSE filter coefficients bf and were derived assuming c„ = c„ for a given set of 

Doppler frequency and delay spread. It is shown that this estimator is robust regardless of 

frequency or time mismatches. With a low Doppler frequency, it has been shown, that a 5- 

tap ( Mf =5) estimator can successfrilly predict the channel. 
[054] To obtain accurate initial channel estimation, a training OFDM block is sent at the beginning 

of a transmission, in which Cj is known to obtain j . The channel parameters for the new 

time epochs and the unknown code words can be successively solved in the forward direction 

(time advance). 

[055] The assumption c„ = c„ caimot be always guaranteed; an incorrectly detected code word 
introduces wrong channel estimation and hence, can cause a wrong detection of the 
successive code words. This error propagation is the dominant impairment of the link 
performance at high SNR. In order to alleviate the problem of error propagation, a training 
OFDM block is inserted every block. In the simulations presented herein, A^^ = 70 is 
considered. 

[056] In FIG. 2, the performance of the coherent reception with RS code and convolutional codes 
for 40-Hz maximum Doppler frequency is shown. The convolutional codes are shown with 
different constraint length (K) ranging from 3 to 9 as dashed lines. The performance of 
convolutional codes is substantially better than that of the RS code, which is shown as a solid 
line with cross lines. In order to achieve a Word Error Rate (WER) of 10"^, the K=9 CC 
needs 4 dB lower SNR than the RS code. In the use of WER in the present invention, it is 
assumed that a word is a codeword. Moreover, the performance of the K=9 CC with channel 
estimation is very close to the one with the ideal channel information shown as a solid line. 

[057] In FIG. 3, the performance at 200-Hz Doppler is shown. In comparison with the RS code, 
the CC's are still superior although the degradation with respect to the idealized case is 
higher due to poorer channel tracking. In fact, an error floor at the high SNR region exists 
due to tracking errors. Once again the RS coded signal is indicated as a solid line with cross 
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lines. The convolutionally coded signals are indicated by dashed lines and ideal channel 
information is indicated by a solid line. 

[058] In FIG. 4, the performance of the K=9 CC with different maximum Doppler frequencies is 
shown. With a low Doppler frequency, a 5-tap ( =5) estimator used here can successfully 
predict the channel and the performance is very close to that with idealized channel 
estimation. However, when the fading is relatively fast, it is difficult to estimate the channel 
correctly and the WER floors on the order of 10*^ can be clearly observed for at a maximum 
Doppler frequency of 200 Hz. That is, it was found that the original method works well in 
slow fading but degrades significantly in fast fading. Once again the ideal Doppler 
frequency is indicated by a solid line. A Doppler frequency of 200 Hz is indicated by a 
dashed line. A Doppler frequency of 175 Hz is indicated by a dashed line with small circles. 
A Doppler frequency of 150 Hz is indicated by a dashed line with triangles and a Doppler 
frequency of 125 Hz is indicated by a dashed line with cross lines. 

[059] FIG. IB depicts a block diagram for a baseband receiver. It is simplified to the extent that 
only channel estimator and decoder portions of the receiver are depicted. The decoder 
further comprises a plurality of demodulators, as well as a maximum ratio combiner, an 
optional deinterleaver (to match an optional interleaver in the transmitter) and a Viterbi 
decoder, all depicted in FIG. 1 A. However, FIG. IB depicts the structure and corresponding 
connections of the receiver depicted in FIG. 1 A. Note, all x, c , H, H are complex. 
Although it may seem strange to have a complex codeword, c , it is very natural to treat 
coding and modulation as a whole like spatial-temporal coding as discussed in "Space-time 
codes for high data rate wireless communication: performance criterion and code 
construction", by V. Tarokh, N. Seshadri, and A. Calderbank, published in IEEE Trans. 
Info. Theory, vol. 44, no. 2, pp. 744-765, Mar. 1998 or coded modulation "Channel coding 
with multilevel/phase signals," by G. Ungerboeck, published in IEEE Trans. Info. Theory, 
vol. IT-28,no. l,pp. 55-67, Jan. 1982. 

[060] At a time instant n, the channel estimator unit has two tasks. One is to produce the 

channel estimates of current time instant, H/„ , with its input X|,n, X2,n, . . . and 

feedback ^y^.M,Mi.„.M, .H^,,., ...by equation (11). The other task is to produce 
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references Hj „ „ v- for estimate processing in the next time instant by equation (10) 
when c„ is available. At a time instant n, when Hj „ ,63 „ are available, the decoder unit 
produces decoded information c„ by (9). 
[0611 The flowchart of this method is shovm in FIG. IC. The related channel estimation 

method is first initialized at step 165-1. Transmitted signals are received at step 165-2. A 
determination is then made at step 165-3 as to whether the received block is a training block. 
If the received block is a training block then c„ is known and 

H„,„ =argmin^||x, „ -H, „c„ ||' 
is calculated at step 165-5. This is a reference for the channel estimation 

/=i 

which is calculated next at step 165-6. The block number is incremented at step 165-7 and a 
determination is made if the end of the frame has been reached at step 165-8. If the current 
block is not a training block then Cn is decoded at step 165-4 
c„ = argmin^lK ^ -H„^c„ ||^ 

m 

is calculated before calculating the reference and channel estimation. 

[062] A method and system using two channel estimation schemes, which improve channel- 

tracking capability based on simplification of the ideal or optimal detector are now presented. 

[063] The present invention uses a COFDM system with convolutional coding (CC). 

Moreover, the present invention uses two iterative processing techniques using current and 
future tentative data decisions in the backward direction for channel estimation and the final 
decoding. This permits joint chaimel estimation and data decoding with improved channel- 
tracking capability, resulting in reliable link performance even under high user mobility 
and/or high RF carrier frequency. System level performance such as retransmission 
probability and packet delay in a system with aggressive fi-equency reuse using dynamic 
packet assignment (DP A) is also presented. 

11 



Chuang 2000-0071 

[064] Since the wireless channel is correlated in time, the optimal joint detection requires 

processing the received signals and the decoded data in the past, current and future epochs. 
From (10) and (1 1), it can be seen that the joint detector in the related is sub-optimal because 
it attempts to predict the channel response at a time instant by only using received signals and 
decoded information in the past. In order to improve channel estimation by applying the 
current information, e.g., decoded data and newly received signal, an iterative solution can be 
applied. 

[065] At each time epoch n , after performing (9), (10) and (11), the present invention does 

not proceed to the time epoch n + 1 immediately. Instead, the present invention uses the 
newly generated , which is derived based on the current information, to substitute for 

the current channel estimate „ and the computations in (9), (10) and (1 1) are repeated. In 

other words, the same procedures are performed twice. This additional iteration allows 
channel estimation to be performed based on the current information. 
[066] In FIG. 5, the improvement of the system performance is shown. At each time instant 

n , if the perfect past chamiel were known, i.e., [H^ ,,,,H^„_,] = [H^ ,,,,H^ in (1 1), the 

new iterative processing can perform within .3 dB from the case in which channel is known. 
On the other hand, even if the perfect past channel information were available, the system 
performance would still be far fi*om the optimal one with an irreducible error floor in the 
related sub-optimal solution. The SNR requirement of the new iterative approach at a WER 
of 10"^ is 1 .2 dB lower than that of the original one. System performance for the original 
signal is indicated by a dashed line. Using an iterative approach is indicated by a dashed line 
with triangles. Using the original approach and having perfect past channel information is 
indicated by a dashed line with circles. Using an iterative approach and having perfect past 
channel information is indicated by a dashed line with cross lines. The ideal estimate is 
indicated by a solid line. 
[067] More than two iterations were attempted but the improvement was limited. 

Therefore, it was determined that another possible improvement probably should come from 
utilizing/w/wre information as described next. 
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Following the insights of the optimal approach, by taking into account both the 
current and the future information another step is performed. To take advantage of the future 
information in the channel estimation process, the FIR filter shown in (11) can be processed 
in the time-reversed fashion, i.e., 

iB,d(H„,.,.,)-d(H„,„.,) = 0. (13) 

Hence, the iterative approach is performed backwards. Here, assume that the iterative 
approach mentioned in the last sub-section has been already carried out from time instant 1 to 
N . Therefore, x„ and c are available with which to perform the iterative approach 
backward at a later time for tentatively detected previous signals; this method is thus referred 
to as iterative backward-processing approach. In order to process backward, H„ „ is stored 
when the decoding resuU is correct (assuming an error detection mechanism is available) and 
c and x„ when a decoding error occurs at a time instant n . 

In general, there are three choices to perform the backward processing. The first 
choice is to process backward after all OFDM-block data are processed. In reality, it is 
not a good choice because the channel memory length K„ is not infinite and the delay of this 
approach may not be acceptable if is large. A feasible approach is to process backward 
after having M, (FIR tap length) consecutive correct code words (OFDM blocks). This 
ensures that the backward processing will start from a reliable channel estimate. The last 
choice is to process backward starting at any other points, which is determined by memory 
and delay requirements. Performance degradation is expected. However, if the backward 
processing is started from a known training OFDM block, lower degradation can be 
achieved. With the last two choices, it is only necessary to store (« N in general) 

OFDM-block data in the memory. 

The iterative processing and iterative backward processing approaches share the 
similar baseband processing block diagram (FIG. 1 A) with the related/original approach. 
The only difference is that the feedback is H,,fr/i^ti,H2^^^,,...,H, „,H2„. and 

H,,„.w,,H2„.^,,...,H,„^,,H2„^, , respectively, when performing iterative processing and 
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iterative backward processing. Their corresponding flowcharts are shown in FIG. 6A and 
FIG. 6B. 

The iterative channel estimation method depicted in FIG. 6A is first initialized at step 
605. The iterative processing for estimating channel characteristics is performed by using the 
system depicted in Figs. 1 A and IB and as described above. Transmitted signals are 
received at step 610. A determination is then made as to whether the received block is a 
training block at step 615 . If the received block is a training block then c„ is known and 
H„ „ = argmin J]||x„ „ - H„ „c„ ||^ is calculated, which is a tentative reference signal, by 

first tentatively decoding the block of the received multicarrier signal at step 620. The 
tentative reference signal is then used to generate a tentative estimation for the channel at 
step 625 given by the equation 

Mi 

XB,d(H„,„_,)-d(H„,,,) = 0. 

i=\ 

The tentative reference signal is then used to generate a tentative estimation for the channel 
at step 625, the block number is incremented at step 630, and a determination is made if the 
end of the frame has been reached at step 635. If the end of the frame has not been reached 
then another block of the received multicarrier signal is accepted for processing at step 610. 
If the current block is not a training block then c„ = argmin^|lx„^ -H„^c„ ||^ is calculated, 

which is a reference signal, by first decoding the block of the received muhicarrier signal at 
step 640, the matrix H„ „ = argminJ^Hx^ „ - H„ „c„ 1|^ is calculated at step 645, and it is then 

used to generate an estimation of channel characteristics given by 
SB/d(H„.„.M)-d(H„.„.,) = 0 

/=1 

at step 650. The block of the received multicarrier signal is then re-decoded using the 
estimation of the channel characteristics at step 655. The method then proceeds to step 620. 
This effectively repeats the reference and channel estimation. The reference and channel 
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estimation are repeated in order to improve the calculations with the tentative reference and 
channel estimation calculations. 
[072] The iterative backward channel estimation method depicted in FIG. 6B is first 

initialized at step 660. The iterative backward processing for estimating channel 
characteristics is performed by and using the system depicted by Figs. 1 A and IB and as 
described above. Transmitted signals are received at step 665. A determination is then made 
as to whether the received block is correct at step 670. If the received block is correct then £„ 
is known and 

H„.„ =argmin5]|lx„ „ -H„ „c„ ||^ 

is calculated, which is a tentative reference signal, by first tentatively decoding the block of 
the received muhicarrier signal at step 675. This tentative reference signal is then used to 
generate a tentative estimation for the channel at step 680 given by 

iB,d(H„,„,,.,)-d(H„.„.,) = 0. 

/=i 

[073] The block number is decremented at step 685 and a determination is made if the 

beginning of the frame has been reached at step 690. If the beginning of the frame has not 
been reached then another block of the received multicarrier signal is accepted for processing 
at step 665. If the current block is not correct block then 6„ 
c„ = argmin2]l|x„.„ -H„^c„ ||^ 

is calculated, which is a reference signal, by first decoding the block of the received 
multicarrier signal at step 692. This reference given by 
H„ „ = argmin2]||x„ „ -H„_„c„ f 

at step 694 is then used to generate an estimation of channel characteristics at step 696 given 
by 

lB,d(H„.„,,.,)-d(H„.„.,) = 0. 



The block of the received multicarrier signal is then re-decoded using the estimation of the 

channel characteristics at step 698, and control is returned to step 675. This effectively 
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repeats the reference and channel estimation. The reference and channel estimation are 
repeated in order to improve the calculations with the tentative reference and channel 
estimation calculations. 

[074] In FIG. 1, the performance of this approach (iterative backward processing) is shown 

for 200 Hz maximum Doppler frequency. In the simulation, the maximum is set to be 
200 and this corresponds to 40 ms. However, due to the low error probability in the high 
SNR region, much shorter storage is required. For instance, the maximum required is 
about 50 at the 5 dB SNR. It is found that nearly optimal performance is achieved with 
iterative backward processing. System performance using the original approach is indicated 
by a dashed line. System performance using the iterative approach is indicated by a dashed 
line with triangles. System performance using the iterative backward processing approach of 
FIG. 6B is indicated by a dashed line with cross lines. The ideal estimate is indicated by a 
solid line. 

[075] In FIG. 8, the performance of this approach with different maximum Doppler 

frequencies is shown. The system still performs well in an environment with maximum 
Doppler frequency as high as 400 Hz. Once again the ideal estimate is indicated by a solid 
line. System performance at 500 Hz using iterative backward processing is indicated by a 
dashed line with crosses (or "x"s). System performance at 400 Hz using iterative backward 
processing is indicated by a dashed line with small circles. System performance at 300 Hz 
using iterative backward processing is indicated by a dashed line with triangles. System 
performance at 200 Hz using iterative backward processing is indicated by a dashed line with 
cross lines. System performance of the original approach at 200 Hz is indicated by a dashed 
line. 

[076] By applying the simple iterative estimation, the impact of noise is reduced by an 

additional round of filtering with newly available information. However, it may not 
eliminate the impact of decoding errors. Using iterative backward-processing, there is a 
chance to re-estimate the chaimel with fewer decoding errors and hence, better performance 
is achieved. 
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[077] It should be noted that the FIR estimator coefficients in the related sub-optimal 

method were optimized for the scenario without any iterative process. It should be noted, 
however, the detector architecture introduced here can be combined with any filter design 
for better tracking. 

[078] To show the robustness of the iterative approaches, a simple averaging 5-tap finite 

impulse response FIR filter (6, = . 2) for the time domain filtering is now considered. As 
shown in FIG. 9, even with these simple FIR coefficients, the method of the present 
invention still outperforms the original method with the FIR that was optimized for a 
particular set of maximum Doppler frequency and delay spread. Therefore, the iterative 
backward-processing approach of the present invention is relatively robust against the 
mismatch between the FIR coefficients and the true channel. The ideal estimate is indicated 
by a solid line. System performance using the original approach and an optimal 5-tap FIR at 
200 Hz is indicate by a dashed and dotted line. System performance using the original 
approach and an optimal; 5-tap FIR at 40 Hz is indicated by a dashed line. System 
performance using iterative backward processing and an optimal 5-tap FIR at 200 Hz is 
indicated by a dashed line with crosses (or "x"s). System performance using iterative 
backward processing and an optimal 5-tap FIR at 40 Hz is indicated by a solid line with 
crosses (or "x"'s). 

[079] Consider a MAC layer configuration to characterize the system-level performance 

under frequency reuse using the improved detection methods. A simulation system of 36 
base stations arranged in a hexagonal pattern is used, each having 3 sectors, with a 20-dB 
ft-ont-to-back ratio and idealized antenna pattem. The same channel can be used everywhere, 
even in different sectors of the same base station, as long as the Symbol to Interference Ratio 
(SIR) in the DPA admission process exceeds 7 dB. The simulation of the simultaneous use 
of the same spectrum by users in different cells results in interference between cells. 
Interference is treated as if it behaves like noise. 

[080] A channel is defined to be a combination of time slot and RF carrier, each consisting 

of 120 subchannels described previously. By way of example, the results for the case with 3 
RF carriers and 8 time slots are now presented. This occupies a total spectrum of less than 



Chuang 2000-0071 

2.5 MHz, including guard bands and other overhead. Each time slot consists of 10 OFDM 
blocks, i.e., 2 msec. One of these OFDM blocks is used for training, as discussed previously, 
while an additional block can be allocated for guard time between time slots. In addition, 
assume that a control slot of duration 4 msec is inserted in the beginning of every frame of 8 
traffic slots to enable paging, assignment and pilot transmission that are required for the DPA 
process. With this conservative assumption of overhead in time and frequency domains, 48 
kb/s (8 OFDM blocks or 960 data bits transmitted in 20 msec) can be delivered using each 
time slots. Once paged, a mobile station (MS) measures the pilot signals to determine the 
desired traffic slots and reports the list back to the base station (BS). The BS then assigns 
traffic channel(s) and informs the MS this assignment for traffic packet delivery. 

[081] Based on the downlink frame structure, four adjacent BS's form a reuse group and 

they take ttims performing the DPA procedure once every 4 frames. The time-reuse groups 
in the entire service area are pre-planned in a fixed and repeated pattern. 

[082] For the propagation model, the average received power decreases with distance d as 

d"^ and the large-scale shadow-fading distribution is log-normal with a standard deviation of 
10 dB. A data-service traffic model, based on wide-area network traffic statistics, which 
exhibit a "self-similar" property when aggregating multiple sources, was used to generate 
packets. 

[083] Automatic Retransmission reQuest (ARQ) is employed for retransmission when a 

packet is received in error. A packet in this case is theoretically 8 code words in each time 
slot, but the error probability is represented by using the WER curves. Since the error 
probability of the 8 code words in a time slot is highly correlated and additional coding is 
usually included for the entire packet, this approximation provides reasonable performance 
estimation for the MAC layer. If a packet cannot be successftiUy delivered in 3 seconds, 
which may be a resuU of traffic overload or excessive interference, it is dropped from the 
queue. The control messages are assumed to be error-free in the designated control slots. 

[084] First consider K=5 and 40 Hz or K=9 and 125 Hz based on the related/original (non- 

iterative) method, both give similar link performance (see Figs. 2 and 5), for a comparison 
between RS and CC coding methods. FIG. 10 shows the average probability of packet 
retransmission, as a function of occupancy for all available (24) channels in each sector. 



Chuang 2000-0071 

This is a measure of QoS (quality of service) experienced by individual users. With a 3-6 % 
target retransmission probability, 1 5-50 % occupancy per radio in each sector is possible with 
this DP A scheme, depending on the use of coding schemes. Clearly, joint channel estimation 
and maximum likelihood detection of CC indicated by a dashed line introduced previously 
provides significant improvement over the case of RS codes, indicated by a solid line with 
cross lines which is also similar to the case of differential demodulation of the RS codes with 
4 transmit antennas. Both results are significantly superior to the efficiency provided by 
current cellular systems, which are typically designed for voice communications with very 
conservative frequency reuse, about 4-7% spectrum occupancy in each sector. Data 
applications, permitting some retransmission delay, and improved link design, introduced 
here, allow much more aggressive and efficiency frequency reuse. 
[085] FIG. 1 1 shows that 1-1 .5 Mb/s can be successfully delivered by each base station with 

an average delay on the order of 40-120 msec. This is a measure for system capacity. It 
indicates that OFDM link and DP A MAC combined enable a spectrally efficient (40%-60% 
b/s/Hz with a conservative assumption of overhead requirements) air interface for broadband 
services, even for the macrocellular environment considered here. Adaptive modulation has 
not been considered in this study, and its use is expected to improve efficiency beyond 1 
b/s/Hz per base station even imder aggressive frequency reuse. The OFDM technology 
discussed herein can provide robust performance with peak-rates scalable with the available 
bandwidth. RS codes are indicated by a solid line with cross lines. CC codes are indicated 
by a dashed line. 

[086] Next, consider the case of high maximum Doppler frequency (200 kHz) and K=9 

(WER curves in FIG. 6) for comparison between sub-optimal ("original") detection method 

discussed earlier and near-optimal iterative backward-processing method presented. FIG. 12 

shows that retransmission probability using the improved method can work well even under 

high maximum Doppler frequency. As a result, QoS can be improved even for high mobility 

users or when higher canier frequency is employed. System performance using the original 

detection method is indicated by a solid line with cross lines. System performance using the 

iterative backward processing approach is indicated by a dashed line. On the other hand, the 

capacity difference is relatively smaller, as shovm by the delay-throughput curves in FIG. 13. 

19 



Chuang 2000-0071 

System performance using the original approach is indicated by a solid line with cross lines 
and system performance using the iterative backward processing is indicated by a dashed line 
This is because both methods give very good radio link performance and the delay is 
dominated by sharing limited number of traffic slots, which is independent of the WER 
performance. In this case, better traffic resource management, such as improved admission 
control, could achieve capacity improvement. If higher Doppler frequency, e.g., 400 Hz, 
were encountered, the improved link performance introduced by the iterative backward- 
processing method, as shown in FIG. 8, would also result in system capacity enhancement. 
[087] Returning to FIG. 1 A to show the detail of the decoder (DEC) unit and its relation 

with the Channel Estimator (CE) unit, with channel estimates „ , c„ can be calculated by 
maximum likelihood (ML) decoding (9), 

c„ =argmin2^K „-H^^c„ ||^ 

[088] In FIG. 1 A, this decoding process is divided into 4 parts, e.g. QPSK Demodulator, a 

maximum ratio combiner (MRC), Deinterleaver and Viterbi decoder. How to separate these 
4 units is now explained. 

[089] We have 

IK„-H„,„cjr 

m 

k m 

= Y.(L\Kn, IP -2Re{x:,,„.,/i„,„.,c„,,}+||/z„,,,c„., n 

k m 

Note that „ ||^ makes no contribution to our minimization and due to binary 

k m 

convolutional code with QPSK modulation, || c„ ^ 1|^ is a constant. ^^SH ^'».'>.* "^'I ^"M H^) 

k m 

makes no contribution to our minimization either. So we concentrate on 



Re{c„,,}- 
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The demodulator and the MRC calculate Y,^^{x'mMK.n.k} and Y,lm{xl,,„,kK,n^) • And 

the remainder is done by the deinterleaver and the decoder. 
[090] Conceptually, the MRC does only energy combining. So, we perform the following 

calculations: 

SRe(^l.">'".".*> = S^^^^^"."^ arg(/7„.„.,)} II II 

m 

m ^ 

Then, Re{x:.„.,arg(;i„.„,,)} and Im{x:^,, arg(A„,„,,)} are actually QPSK demodulation. And 

the remainder is the MRC. 

[0911 The radio-link performance of our COFDM system with the novel and non-obvious 

combination of two channel estimation schemes under fast fading validate the design 
methodology and system for near optimal joint channel estimation and data detection. 

[092] For a COFDM system, coherent detection can significantly improve radio link 

performance, but it requires a channel estimator to achieve the potential gain. The present 
invention applies a forward channel estimation and maximum likelihood decoding method 
for a COFDM system with convolutional codes and a simple interleaving scheme. It is 
shown that the performance of such a system significantly outperforms the one with Reed- 
Solomon codes. The performance under different fading rates was studied to understand the 
performance limits and areas for improvement. Simulations showed that this detection 
method is able to achieve near optimal performance in a wireless environment with 
maximum Doppler frequency as high as 100 Hz. However, as fading rate increases fiirther, 
irreducible error floor is introduced by the limitation in channel-tracking capability. A near- 
optimal detection method has been presented to improve the channel tracking performance, 
which is based on iterative processing and iterative backward processing of channel 
estimation and data decoding. Within practical values of interest, error floor can be 
eliminated even for very high fading rates, which is important for cases with high mobility or 
when higher carrier frequencies are considered in the future. Radio system performance 
based on the improved link techniques and a medium access control protocol using dynamic 

packet assignment was also studied. With radio link improvement, system throughput and 
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delay can be significantly enhanced for better quality of service and spectrum efficiency. 
The system is also greatly simplified by allowing reuse factor of one without frequency 
planning. 

It should be clear from the foregoing that the objectives of the invention have been 
met. While particular embodiments of the present invention have been described and 
illustrated, it should be noted that the invention is not limited thereto since modifications may 
be made by persons skilled in the art. The present application contemplates any and all 
modifications that fall within the spirit and scope of the underlying invention disclosed and 
claimed herein. 
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